Superconducting strip in an oblique magnetic field 
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As an example for a seemingly simple but actually intricate problem, we study the Bean critical 
state in a superconducting strip of finite thickness d and width 2w S> d placed in an oblique magnetic 
field. The analytical solution is obtained to leading order in the small parameter d/w. The critical 
state depends on how the applied magnetic field is switched on, e.g., at constant tilt angle, or first 
the perpendicular and then the parallel field component. For these two basic scenarios we obtain 
the distributions of current density and magnetic field in the critical states. In particular, we find 
the shapes of the flux-free core and of the lines separating regions with opposite direction of the 
critical currents, the detailed magnetic fleld lines (along the vortex lines), and both components 
of the magnetic moment. The component of the magnetic moment parallel to the strip plane is a 
nonmonotonic function of the applied magnetic field. 

PACS numbers: 74.25.Qt, 74.25.Sv 
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I. INTRODUCTION 

Platelet-like type-II superconductors in a magnetic 
field applied at some angle to the normal of their plane 
are frequently investigated in various experiments, see, 
e.g., Refs. |l||2y|jl5l However, even for the simplest case 
of an infinitely long strip placed in an oblique magnetic 
field, the critical state was theoretically studied only in 
the situation when the magnitude of the applied mag- 
netic field Ha considerably exceeds the field of fuU-fiux 
penetration into the sample^ Hp. The attempt to inves- 
tigate the critical state in fields Ha < Hp led to incorrect 
results)^ since an essential feature of this state was over- 
looked, as will be evident from our analysis in Sec. II 
B. 

In this paper we consider the following basic situa- 
tion: A thin superconducting strip fills the space \x\ < w, 
\y\ < oo, \z\ < d/2 with d <C w; a constant and homoge- 
neous external magnetic field Ha is applied at an angle 9 
to the z axis (^Hax — Ha sin^, Hay — 0, Haz — Ha cos^). 
It is assumed that the thickness of the strip, d, exceeds 
the London penetration depth, the critical current den- 
sity jc does not depend on the local induction B (Bean 
modetfiS), and the lower critical field Hd is sufficiently 
small so that we may take B = hqH. We consider 
two scenarios of switching on the external magnetic field: 
First, the magnitude of the external field increases from 
to Ha at a fixed angle 9; second, one turns on Haz 
first and then Hax- Interestingly, these scenarios lead to 
different critical states. 

Takinginto account the result of Ref. (see also 
Refs. llOlllll) . the smallness of the parameter d/w enables 
us to split the two-dimensional critical state problem 
for the strip of finite thickness into two simpler prob- 
lems: A one-dimensional problem across the thickness of 
the sample, and a problem for the infinitely thin strip. 
This splitting becomes possible since under the condition 



d/w 1 the magnetic fields and currents in the criti- 
cal state essentially change along the x direction only on 
scales which considerably exceed the thickness d. 

The solution of the critical state problem for the in- 
finitely thin strip is knowniiSiiii^ The z component of 
the magnetic field is completely screened by the currents 
fiowing in the region |a;| < a, i.e., one has Hz = there. 
The length a is described by the simple formula: 



a{h) 



cosh(ft,cos6') 



(1) 



Here and below h = Ha/ He, He — Jc/t^-, Jc = jcd, and w 
is taken as the unit of length {w — 1). In this region of the 

strip, \x\ < a, the sheet current J{x) — ^'^^^^r^j{x, z)dz 

(with current density j along y) is given by 



2 ^ x^/l 
J[x) = Jcarctan- 



(2) 



On the other hand, at a < |x| < 1 where Hz{x) ^ 0, one 
has J = — sign(x)Jc with sign(a;) — ±1 for x > and 
a; < 0, respectively. The explicit form of Hz{x) in this 
region of the strip is presented in Appendix A. 

In the region |a;| < a of the real strip (with d ^ 0), the 
fiux lines are practically parallel to the planes of the strip 
and penetrate into the sample across its thickness from 
the upper and lower surfaces of the superconductor. The 
penetrating flux fronts form the boundary of a flux-free 
core, z^{x), which thus consists of an upper and a lower 
branch. Below we consider only the upper branch since 
from symmetry considerations one has: — zi^°™'^(— x) = 
z7P"(x) = z^{x). 

Following our idea of splitting the critical state prob- 
lem, we consider a small section of the strip around an ar- 
bitrary point X {\x\ < a) as an "infinite" slab of thickness 
d placed in a parallel dc magnetic field Hax and carrying 
a sheet current J{x), Eq. The critical state in such 



2 



a slab is well knowrifiii^ and this enables us to find the 
flux fronts, the distribution of the magnetic fields, and 
the currents across the thickness of the strip in the re- 
gion |a;| < a. Since the critical state in the slab depends 
on how Hax and J was turned on, the above-mentioned 
dependence of the critical state in the strip on the pre- 
history of Hax and of Haz appears. Of course, a similar 
procedure may be used in the region a < \x\ < 1 to find 
the distribution of the magnetic fields, but there the ap- 
propriate analysis is trivial since j{x,z) is constant; we 
thus do not discuss it below. 



II. MAGNETIC FIELD IS INCREASED AT 
CONSTANT TILT ANGLE 

In the case of the first scenario of switching on the 
magnetic field when 6 =const and h has increased mono- 
tonically, it is convenient to introduce the function. 



F{x,h) 



h sin 6 — arctan 



(3) 



which is proportional to the x component of the magnetic 
field on the upper surface of the strip at < a, H^^s(x) = 
Hax + 0.5 J {x) — {Jc/2)F{x, h). Below we shall also use 
the two characteristic fields: 



K = 



2 sin 6* ' 

and hf defined by the equations: 

arctan[tan 6 th(u cos 6)] 



cosh(i 



sin 9 cosh(/i / cos 9) , 



(5) 



where u is some parameter. The meaning of these fields 
will become clear below. 



A. Interval < h < hf 

Consider first the flux front in the interval < h < hf. 
It is essential that there exists a point on the upper plane 
of the strip where the derivative dH^six) / dx vanishes. A 
simple calculation gives that this occurs at the point with 
the coordinate xi, 



x\ (h) = aih) sin ( 



(6) 



Thus, when h increases, the flux lines at —a{h) < x < 
xi(h) monotonically penetrate into the strip through its 
upper surface, and the shape of the core in this interval 
of X is determined by the equation: 



z^{x) ^ [dl2)\V- F{x,h)\. 



(7) 



This formula follows from the well-known distribution 
of the magnetic field in the "slab" shown in Fig. 1 



("Bean profiles"). On the other hand, in the interval 
x\(K) < X < xo{h) the flux lines leave the sample, and at 
X > xq vortices of opposite sign penetrate into the strip. 
Fig. 1. Here the point xo{h) is found from the condition 
Husixo) = 0, i.e., from 



FixQ,h)=0. 



(8) 



It is clear from the inspection of Fig. 1 that in the interval 
xi{h) < X < X2Q1) the shape of the core is determined 
by the flux front occurring at the field /i* which has to 
be found from the equation 



X = a{h^) sm9. 
Thus, in this interval one has 

z^{x) = {d/2)[l~F{x,h,)]. 



(9) 



(10) 



In the same interval there is also a front zi{x) separating 
the regions of the strip with opposite signs of the critical 
current density, see Figs. 1,|21 This front is described by 
the formula: 



zi{x) = (d/4)[2 + F{x, h) - F{x, h,)] . 



(11) 



At the point X2{h) the front zi{x) reaches the boundary 
of the core Zj(x), and hence, this X2{h) is determined by 
the condition z^{x2) — zi{x2). Using formulas l|^- Hll|l . 



(4) one then obtains for X2'. 



X2 — a(u) sin 9, 
F{x2,h) + F{x2,u)^0. 



(12) 



At X2{h) < X < a{h) the critical current density has only 
negative sign, j — —jc, at z > z^, and we arrive at 



z^{x) = {d/2)[l + F{x,h)]. 



(13) 



We see that in the interval of the magnetic fields < 
h < hf the width of the flux- free core is equal to 2a{h), 
but its size along z is less than d, see Fig. 2. When 
h ^ hf, the difference a{h) — X2{h) tends to zero, and at 
h = hf one has X2{h) = a{h). With the use of Eqs. ^ 
and \i2\ this condition may be rewritten as Eqs. ((S)). 



B. Interval hf < h < hp 



When hf < h < hp, the upper and the lower branches 
of the flux-free core merge in the intervals X3 < \x\ < a, 
and hence the size of the flux-free core in the x direction, 
2x3{h), becomes less than 2a{h), see Figs. 3, 4. Here 
X3{h), which lies between xi{h) and a{h), is determined 



by the condition, z^^PP'^'^ixa) 
plicit form by the equations: 



ylowcr 



(X3), or in the ex- 



^3 — a{u) sin0, 
F{-X3,h) + F{x3,u) = 2, 



(14) 
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FIG. 1: Bean profiles of the magnetic field Hx{z) across 
the strip at 5 positions x in the strip at h < hf. Here 9 = 
45°, Ha = 0.6Hp < Hf ^ 0.834Hp, thus xi = 0.478, xo = 
0.597, X2 = 0.617, a = 0.677 in units w (see Fig. 6). The 
characteristic z values z.y{x), zi{x) = —Zry{—x), and zi{x) are 
defined in the text and in Fig. 2. Shown are: a) < a; = 
0.35 < xi, h) xi < X — 0.57 < xq, c) a;o < 2; = 0.61 < X2, d) 
a;2 < 2; = 0.65 < a, e) a; > a = 0.677. For all these profiles, 
2Hx/Jc is F{x,h) on the upper surface, and F{—x,h) on the 
lower surface. The continuation of the increasing parts of the 
profiles b) and c) intersects the upper surface at F{x,h,). 



where we have taken into account the symmetry of the 
flux-free core and formulas 0, 1(1171) . In other words, 
we find that not only the z-size of the core is less than 
d, but also its a;-size is less than the width of the re- 
gion in the strip, 2a{h), where Hz = 0. In the inter- 
val —x^^h) < X < xi{h), the core is again described by 
Eq. Q, while in the interval xi{h) < x < X3{h) it is given 
by Eqs. lO, (|10|l . When h becomes equal to hp defined 
by formula Q), the point X3 reaches xi, and moreover, 
the core disappears at all since at this field the difference 
z^^PP'^'ix) - z\°'^'''{x) vanishes even for |a;| < xi. Thus, 
hp — 0.57r/sin6' is the field of full penetration of flux 
into the strip in the oblique magnetic field. This field 
has a simple meaning. In usual units one finds for the 
X component of the penetration field, Hchp sin 9 = Jc/"^, 
i.e., the penetration occurs when the x component of the 
applied field has completely penetrated into the slab. 

Interestingly, the part of the boundary of the core in 
the interval xi{h) < x < X3{h), as well as its part in the 
interval xi{h) < x < X2Q1) for the fields < /i < /i/, is 
described by a universal function of z on a; which does 
not depend on h at all, see Eqs. J^, (|10|) . The upper 
corner of the core, i.e., the point X2{h) ior < h < hf 
or the point X3{h) for /i/ < h < hp, moves just along the 
line described by this function when h increases. Fig. 5. 

As to the front separating the regions of the strip with 
opposite signs of the critical current density, it is still 
described by Eq. Hll|l in the region xi < x < a even for 
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FIG. 2: The flux and current fronts in a thin strip in an 
increasing magnetic field Ha that is inclined by a constant 
angle 9 — 45° away from the normal of the strip plane. Shown 
is the state when Ha = 0.7Hp < Hf = Q.834Hp where Hp = 
(Jc/2)/sin0 is the field of full penetration, see text. The 
current-free core is delimited by the functions z-/{x) (bold 
line) and —z.y{—x). The current density to the left of this 
core and of the "tails" z-i_{x), —zi{—x), is j = jc, and to the 
right j — —jc- The lower plot enlarges the region near the 
points X = xi, X — X2, and a; = a on these curves, see text. 
The two dashed lines are extensions of the fronts depicted 
in the regions —a < x < xi and X2 < x < a and cut the 
upper surface 2: = d/2 at the point x — xq where the local Hx 
vanishes, Eq. (8). 

hf < h < hp. However, in the interval x^ < x < a, apart 
from this upper branch of the front, zi{x), a lower branch 
Z2{x) appears, and these branches join each other with 
vertical slope at x = a{h), Fig. 4. Knowing the upper 
branch, one can find the lower branch from the given 
value of the sheet current J{x), Eq. JSjl, yielding 

Z2ix) = {d/4)[F{-x, h) - F{x, K)] , (15) 

where h^(x) is found from 

X = a{ht,) smO. 

A complete set of flux and current fronts is shown in 
Fig.inifor three tilt angles 6* = 30°, 45°, and 60°. Note 
that when Ha increases, not only does the flux-free core 
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FIG. 3: Bean profiles of the magnetic field Hx{z) across 
the same strip as in Fig. f but for higher applied field h > 
hf, Ha = 0.9SHp > Hf = 0.834J/"p, thus xi = 0.3II, = 
0.422, xo = 0.438, a = 0.440 in units w (see Fig. 4). The 
characteristic z values z^{x), zi{x) = —Zj{—x), zi{x), and 
Z2{x) are defined in the text and in Fig. 4. Shown are the 5 
profiles: a) < a; = 0.15 < xi, h) xi < x = 0.40 < x^, c) 
X = X3 = 0.422, d) X3 < X = 0.43 < a, e) x > a ^ 0.440. 
Profile d) does not cut the field- free core and has Hx(z) > 
everywhere. The values of Hx on the upper and lower surfaces 
are given by the same expressions F{x, h) and F{—x, h) as in 
Fig. 1. 



shrink but also the current front separating the regions 
with j = ±jc shifts in the sample. In other words, the 
current distribution changes not only near the core but 
also in the region away from it. This result disproves the 
main assumption of Ref. that the currents can change 
only at the flux front but the critical currents remain 
unchanged in the regions penetrated by flux lines. In 
Ref. the flux fronts in inclined fleld are thus incorrect. 



C. Region h > hp 

Although at h > hp the z component of the magnetic 
field does not penetrate into the region |a;| < a{h) of 
the strip, its x component completely penetrates into the 
sample, and the flux-free core is absent. Fig. 7. The two 
branches of the boundary separating the regions of the 
strip with opposite signs of the critical current density 
are still described by Eqs. (|11|) and IjlSfl in the interval 
xi < X < a. In the interval |a;| < xi only one of these 
branches exists, which is the continuation of the Z2{x) 
and is described by the formula: 



Z2ix) 
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FIG. 4: The fronts as in Fig. 2 but at a larger field Ha = 
0.93-ffp > Hj. As in Fig. 2 the flux-free core is formed by the 
curves z~f(x), —z~^{—x\ but the tail separating regions with 
j = ±jc is now composed of two functions 2:1(2:) and 22 (x), 
which join vertically ai x — a and reach the core at x — x-^. 



This formula follows from the given value of the sheet 
current J{x), Eq. Q. 

Interestingly, even at high magnetic fields h ^ hp the 
current front is S-shaped and does not tend to a straight 
line as it was assumed in Ref. 0. Probably, it is for this 
reason that there is a disagreement between the theoret- 
ical and experimental results in Ref. at 9 ~ 7r/2. How- 
ever, it is necessary to keep in mind the following: Our 
approximation based on the splitting procedure is valid 
if dz^/dx <ti 1 and dzi^2/dx <C 1. These inequalities are 
fulfilled almost everywhere in the strip when the charac- 
teristic scales in the x direction (i.e., si, and a — xi) con- 
siderably exceed the thickness d. Thus, the region h > hp 
may be considered within our approximation only when 
(7r/2)cot6' < ln(2ew/c?), i.e., when the angle 9 is not too 
small, see inset in Fig. 7. Otherwise, the component Hz 
completely penetrates into the sample at lower Ha than 
the component H^ does. In this case the field of full 
penetration is Hp/cos9 where H^ = ( Jc/tt) ln(2ew/d) 
is the penetration field at 6* = 0,^^ and the current fronts 
will differ from those shown in Fig. 7. Note that the 
angular dependence of the true penetration field is given 
hy Hp (9) = mm{Hp / cos 9 , Jc/2sm9) and is a nonmono- 
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FIG. 5: Two flux-free cores in a thin strip in increasing mag- 
netic field Ha inclined hy 6 = 45° as in Figs. 2, 4. The bold 
dotted line shows the universal curve, Eqs. (9,10), on which 
the upper part xi < x < X2 ot xs of all flux fronts lies when 
Ha < Hp, see text. The depicted cores belong to Ha/ Hp = 0.5 
and Ha/ Hp = 0.95. The circles mark the characteristic points 
xi, X2, and X3 on the fronts. 



tonic function, see inset in Fig. 7. 



III. MAGNETIC FIELD COMPONENTS ARE 
SWITCHED ON SUCCESSIVELY 

Consider now the scenario when the field Haz is 
switched on first and then one switches on Hax- The re- 
sulting magnetic field again is Hax = Ha sin 6*, Hay = 0, 
Haz — Ha cos 6. In this case the analysis of the critical 
state is quite similar to that presented in the Sec. II, and 
we present only the results here. 

The field of full penetration of flux is still described by 
formula l@J. However, at any h < hp the x size of the flux- 
free core, 2ai, is less than the width 2a = 2/ cosh(/i cos 9) 
of the region where Hz — 0. This ai is determined by 
the formula J(ai) = Hax — Jc, or explicitly, by 



ai 



acos(/ia;/2) 



[l-a2 sin2(/^,/2)] 



1/2 ' 



(17) 



where hx — hsinO — Hax /He- One more characteristic 
scale is xi{h) determined by the relation: J(ii) = —Hax, 
which leads to the explicit expression for xi: 



Xi 



a sin(/ia;/2) 



[l-a'cos^hx/2)f^' 



(18) 



In the interval — ai < x < xi the shape of the flux- 
free core is described by formula (0. But in the region 
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FIG. 6: Flux and current fronts in a thin strip in an increasing 
magnetic field Ha inclined by three angles 6 = 30°, 9 = 45°, 
and = 60° away from the strip normal. Shown are the 
complete fronts for six fields Ha/ Hp = 0.2, 0.4, 0.6, 0.8, 0.93, 
and 1 (bold line), where Hp — (Jc/2)/ sin 6', see also Figs. 1-5. 



Xi < X < ai one has 
z^{x) = ^ 



, 2 xVl ~ 

1 arrta.n 

TT \/0?- — X^ 



(19) 



In this interval of x there is also a boundary z-^ix) sepa- 
rating the regions of the strip with opposite directions of 
jc, see Fig. 8. This horizontal line is described by 



^3(2;) = - 



hx 

1 

TT 



(20) 



This boundary continues in the region ai < a; < a where 
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FIG. 7: The current fronts separating the regions with crit- 
ical currents of opposite direction in a thin strip in increas- 
ing magnetic field Ha inclined hy 6 = 45°. At the depicted 
fields Ha/ Hp — 1, 1.4, 2, 3 above the penetration field Hp — 
(Jc/2)/sin&, the flux-free core has collapsed into one line, 
and the front is composed of the three curves 2:2(0 < a; < xi) 
[Eq. (16)], 2:2(2:1 < X < a) [lower branch, Eq. (15)], and zi{a > 
X > xi) [upper branch, Eq. (11)]. Note that even at large Ha 
these front lines are not straight lines as it might be expected. 
At a <C w these fronts collapse into one curve when plot- 
ted versus x/a. Inset: the dashed line h ^ ln(2ew/d)/ cos 6* 
schematically shows the boundary above which our splitting 
procedure fails, see Sec. II C. The solid line gives the penetra- 
tion field hp{6) = TvHp/Jc = min[ln(2ei(;/ci)/cos0, O.dir/ sin 9] 
and the dash-dotted line shows hf{9), Eq. 



FIG. 8: Current and flux fronts in a thin strip to which 
first the perpendicular magnetic field Haz is applied and then 
the in-plane component Hax is increased from zero to Hax = 
Haz = 0.7Jc/2, resulting in a final tilt angle 9 = 45° (scenario 
2, Sec. III). The solid line from x — —ai via xi, ai, —xi 
to —ai forms a flux-free core. This core is connected to the 
surfaces by tails composed of a horizontal part Z3{xi < x < 
ai), and a curved part 23(01 < x < a), ending at the point 
{a,d/2). The dashed lines give the extensions of the fronts 
which exactly coincide with the dashed lines shown in Fig. 2 
and end at xq. The dash-dotted lines are the current fronts 
of this scenario above full penetration, at the same fields as 
in Fig. 7 for 9 = 45°. Note that these fronts are monotonia, 
while the corresponding fronts of scenario 1 in Fig. 7 are S- 
shaped. 



it is given by an expression coinciding with Eq. (|16|l : 



23(2;) = - arctan^=== . (21) 

TT Va — a; 

When h approaches hp, the point ai tends to xi, and the 
difference z"PP™(a;) — z^™''''{x) vanishes simultaneously 
for all \x\ < xi. Ai h — hp the flux-free core disap- 
pears, while the boundary 2:3(0;) exists at h > hp, and 
it is described by the expression (|21|l in the whole inter- 
val —a<x<a, see Figs. 8, 9. When Hax > >^c/2 is 
increased further, this saturated current front does not 
change any more. 

Thus, we see that the shape of the flux- free core and 
the boundary between the regions with opposite direc- 
tions of the critical current density do not coincide with 
those described in Sec. II and thus depend on the mag- 
netic history. 

IV. MAGNETIC FIELD LINES 

Using the obtained results, it is easy to find the distri- 
bution of the magnetic fields in the critical state of the 



strip. One may either integrate over the current-carrying 
area, noting that each current path has the magnetic field 
of a straight wire (Fig. 10). Or one may derive analytical 
expressions using our splitting approximation (Fig. 11). 
The z component, Hz{x), in this approximation does not 
depend on z and is given by the formulas of R,efs. ll2liai3 
[or by Eq. (A8) in compact form], and the x component 
is 

Hx{x,z)^Hx{x,~d/2)+ [ jy{x,z')dz', (22) 

J-d/2 

where Hx{x, —d/2) is the field on the lower surface of the 
strip. At |a;| < a one has i/^; (a;, -d/2) = {Jc/2)F{-x,h), 
while Hx{x,—d/2) = HaSinO + ( Jc/2)sign(a;) at a < 
|a;| < 1. Taking into account that jy{x,z) has only the 
values jc, —jcj or 0, and knowing Zj{x) and the bound- 
aries between the regions with ±jc, one can easily calcu- 
late Hx{x, z) everywhere in the strip and near the strip 
explicitly. 

As an example, Figs. 10 and 11 show the magnetic field 
lines (parallel to the Abrikosov vortex lines) in the strip 
and near the strip for the first scenario of switching on 
the magnetic field at constant tilt angle 9 till Ha — 0.7 Hp 
is reached. Both figures show the field lines obtained as 
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FIG. 9: The flux and current fronts in a thin strip in which 
first Haz ~ Jc/4: is applied and then Hax is increased from 
zero to full penetration, 8Hax/Jc ~ 0, 1, 2, 3, 4 (scenario 
2, Sec. III). The complete front for Hax = Jc/i is depicted 
as bold line, see also Fig. 8. The front 4 through the strip 
center (0,0) applies to all Hax > >/c/2 and depends only on 
Haz- Note that the depicted fronts correspond to different 
tilt angles 0, e.g., the saturated front corresponds to 6^ = 
arctan(_ffaa;/J^az) > arctan2 « 63°. 



contour lines of the vector potential Ay{x,z) related to 
ll{x,z) = V X (yAy). Figure 10 depicts the field hues 
calculated directly from Ampere's law using the currents 
obtained in Sec. II A and formula (|Aip of Appendix A 
a.t d/w = 0.08. It is important that this current distri- 
bution indeed leads to a flux-free core which is close to 
that obtained in Sec. II A. On the other hand, figure 11 
uses expressions (jA2|l - (|A8|) of Appendix A for the same 
d/w = 0.08. These expressions were derived with our 
splitting procedure. It can be seen that the agreement of 
both field-line patterns is good, but the fine details near 
the current fronts can be more easily resolved in Fig. 11 
(top) which is based on the simple analytical formulas. 
In particular one can see that the field lines exactly fiow 
around the core in which j = 0, and some field lines cut 
the line ( "tail" ) that separates regions with jy = ±jc and 
runs from x = xi on the upper surface to the cusp of the 
core at a; = a;2, see also Fig. 2. The slight wiggle of the 
field lines occurring near = a in Fig. 11 (bottom) is 
an artifact, since the condition for the splitting proce- 
dure, dz^/dx ^ 1, fails at this point. However, a more 
detailed analysis shows that the difference between the 
field-line patterns of Figs. 10 and 11 manifests itself only 
in narrow intervals near x = ±a. 




FIG. 10: The magnetic field lines for a strip of width 2w and 
thickness d = O.OStu in an applied field inclined hy 9 = 45° and 
increased from zero to Ha = 0.7-ffp as in Fig. 2 (scenario 1). 
These field lines were computed as equidistant contour lines 
of the vector potential Ay{x,z), Eq. (Al), yielding a density 
of lines proportional to the local magnetic field H{x,z). The 
gray area shows the current and field-free core. The dashed 
line is along Ha- The upper plot is to scale, the lower plot is 
stretched along 2 by a factor of 4. 



V. MAGNETIC MOMENT 

In an oblique magnetic field, apart from the z com- 
ponent of the magnetic moment of the strip, Mz, an x 
component appears, and both can be investigated 
experimentally^ The expression for Mz (per unit length 
along y) is knownsi^i 



Mz 

JcW^ 



tanh(ft, cos ( 



(23) 



Knowing z^{x) and the boundaries between the regions 
with opposite signs of the critical current density, given 
in Sees. II and III, one can calculate for the strip 
from the formula: 



rd/2 

dx / zjy{x, z)dz . 

-a J-d/2 



(24) 



Here Mx is given per unit length along y, and we have 
taken into account that only the region of the strip, \x\ < 
a, gives a nonzero contribution to Mx- The saturation 
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FIG. 11: The magnetic field lines for the same strip as in 
Fig. 10 (rf = O.OSu;, 9 = 45°, Ha = Q.THp) from Eqs. lIMl - 
l|A7|l . To reveal the details near the current fronts, the 
field lines here are nonequidistant contour lines of the an- 
alytic vector potential Ay{x,y) of Appendix A, with levels 
Ai, (X !/'^sign(!/), u = 0, ±1, ±2, . . ., yielding more lines at low 
fields. The gray area shows the current and field-free core. 
The dashed line marks the surface of the strip and the dotted 
line shows x = xq, Eq. defined by Hx{xo, d/2) = 0. Note 
that field lines cut the line (from a; = a;i to a; = a;2 marked by 
bold dots) that separates regions with jy = ijc- The lower 
plot shows the same case on difi^erent x and z scales. 



value of Mx, which is achieved in high fields for 6 = tt/2, 
is M^'^* = - Jcdw/2, see Appendix B. 

It is known for the infinitely thin strip^*' that whatever 
its length in y, the ends of the strip (in y) always give the 
same contribution to as that caused by the currents 
jy flowing in the y direction. When Mx is calculated, it 
is necessary to allow for the fact that near the ends of the 
strip the currents may have not only y and z components 
but also an x component, i.e., the problem becomes three 
dimensional. However, using the conservation law for the 
current, divj — 0, one can show that even in this three 
dimensional case the ends of the strip strictly double M^- 
It is for this reason that the factor 1/2 was omitted in 
formula (|^ . 

In Fig. 12 we compare the iJa-dependences of for 
the two scenarios of switching on the magnetic field. 
Note that Mx (Ha) of scenario 2 is always larger than 
Mx^\Ha) of scenario 1, except for the trivial angles 
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FIG. 12: Magnetic-moment component Mx{Ha)/Ml'^^ from 
Eq. 12411 and Appendix B, for scenario 1 (solid lines) and 
scenario 2 (dashed lines) plotted versus the applied magnetic 
field Ha in units of the penetration field Hp — Jc/2sm9 for 
tilt angles 8^15°, 30°, 45°, 60°, 75°, and 85°. Here A'/^' = 
-Jcdw/2. The dot-dashed parabola M^/Mf-^ = 1 - (1 - 
Ha/Hpf applies to £» = 90° {Ha along x), and for 6) = {Ha 
along z) one has Mx — 0. For all other angles the \Mx{Ha)\ 
for scenario 2 is larger than for scenario 1. 



6 = where Mx 0, and 6* = 7r/2 where Mx/Mf^^ = 
1 — (1 — Ha/HpY. All Mx{Ha) exhibit a maximum of 
height M^'^'^/MI^* « 26'/7r occurring at Ha/ Hp « 26'/7r, 
and they have the same slope dMx/dHa = 2Ml^'^ / Hp at 



f 2") 

Ha = 0. The difference M^ ' - M^'' is also maximum 
near Hp, see the dashed curve in Fig. 13. Figure 13 plots 
the Mi^-* of scenario 1 as in Fig. 12 but with both ab- 
scissa and ordinate stretched by a factor tt/ {29) > 1 such 
that the approximate scaling of the AIx-^\Ha) at not too 

large Ha/ Hp is seen. The Mx {Ha) curves of scenario 2 
scale even better. 

The nonmonotonic dependence of Mx on Ha can be 
understood from the following arguments: In the region 
|x| < a the X component of the external magnetic field, 
Hax — HaSinO, leads to an asymmetric distribution of 
the currents over z, see Figs. 1-9. It is this asymmetry 
that generates the Mx component. Thus, Mx can be 
estimated as follows: Mx ~ M|^* • {a/w){Ha/ Hp). The 
factor a/w decreases with Ha, see Eq. and its product 
with the increasing factor Ha /Hp leads to the observed 
nonmonotonic behavior of Mx{Ha). 



r(i) 



VI. CONCLUSIONS 

We solve the critical state problem for a strip of finite 
thickness in an oblique magnetic field. Two scenarios of 
switching on the external magnetic field are considered: 
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FIG. 13: The same magnetization curves as in Fig. 12 for 
scenario 1, plotted as Mx/{Ml''^ x 29 /-iv) versus Ha/ {Hp x 
26 /ti). The dots mark the point where Ha = Hp. The dot- 
dashed curve gives the limit 6 = 90°. The dashed line shows 
the difference of the of scenarios 2 and 1 for 9 = 45° in 
form of 15 x (Afi^' - Mi^'^)/{MT^ x 26/-r). 



magnetic field Ha exceeding the penetration field, the 
current front separating the regions of the strip with ±jc 
generally is not a straight line and can shift in the sample 
when Ha is increased further. 

The above general features of the critical sate are ro- 
bust and hold even if the strip is not very thin or if its 
cross section is not rectangular, i.e., these features are 
common for all critical states in nonsymmetric situations. 
However, the fine details of our results (e.g., the short tail 
from xi to X2 in Fig. 2) require that all characteristic 
lengths in the plane of the strip considerably exceed its 
thickness. If some length of the flux-free core does not 
satisfy this condition, one may expect deviations from 
the presented results in this region of the core. Further- 
more, the temperature should be low enough that flux 
creep does not smear these details, i.e., the creep expo- 
nent n in the current-voltage law E{j) cx (j/jc)" should 
be large. A detailed numerical investigation of the effect 
of flux creep is under way. Note that the detailed shape 
of the flux-free core and of the boundaries between re- 
gions with opposite critical current in principle can be 
investigated via the in-plane component of the magnetic 
moment, while it has little influence on the magnetic field 
on the surface. 



(1) the magnetic field is increased at a constant tilt angle 
9, and (2) the magnetic field components are switched on 
successively. The resulting critical states are different in 
these two cases, even after the flux has fully penetrated 
the strip. 

Another characteristic feature of both states is that, 
below the field of full penetration, the height of the fiux- 
free core is less than the strip thickness, i.e., the core 
does not reach the flat surfaces but is connected to them 
by lines ( "tails" ) that separate regions with opposite di- 
rection of the critical currentsiii Moreover, the width of 
the core may be narrower than the region of the strip 
in which the z component of the magnetic field (i.e., the 
component perpendicular to the plane of the strip) van- 
ishes. 

One more interesting feature of the critical states in 
strips in an oblique magnetic field follows from the data 
of Figs. 6 and 9: When the applied field increases, and 
hence the flux lines further penetrate into the sample, the 
current distribution changes not only near the flux front 
but also away from it, i.e., in the regions where the critical 
state was established before. Note that this feature is 
also seen in figures of Ref. in which the critical state 
of a rotating cylinder was considered in magnetic fields 
perpendicular to its axis. These findings mean that this 
property of the current distributions in the critical states 
is characteristic of the general case when the geometry of 
experiment is not too symmetric, while the usual change 
of currents at the penetrating flux front occurs only in 
special symmetric situations (e.g., when the tilt angle 9 — 
0, 7r/2, or when the cylinder does not rotate). Finally, we 
found the somewhat unexpected result that in an oblique 
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APPENDIX A: VECTOR POTENTIAL 

The magnetic field lines of a strip parallel to y coincide 
with the contour lines of the vector potential Ay{x,z) 
related to the current density jy{x, z) by 'S/'^Ay = —jy or 

Ay{v)^ jd\'3y{v')^^^^^, (Al) 

with r = {x, z). 

For the special case of a thin strip in an oblique applied 
field, we can also find Ay{x, z) using our splitting proce- 
dure. Here, as an example, we give the expressions for 
Ay in the case of scenario 1 at ft- < ft-/ (see Fig. 2). Inside 
the core delimited by the two lines z^(x) and —z^{—x), 
one has H^ — H^ ~ jy = 0, and we may put Ay = 
there. Within the core width (|a::| < a) one obtains: 

Ay{x,z)=- f H^{x,z')dz'. (A2) 

Inserting the H^ from Sec. IV into Eq. (jA2|) . one finds 
explicit formulas for Ay. Equivalently, the Ay inside 
the core width can be calculated directly from the equa- 
tion d'^Ay/dz^ = —jy and the current density of Sec. II. 
Eventually, we obtain the following expressions for Ay = 
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{jcl'^) a,y{x, z), depending on a, xi, X2, z^{x), and 21(2;): 
For —a < x < Xi, < z < d/2 and for xi < x < X2, 

Zj < Z < Z\. 



(A3) 



for Xi < X < X2, zi < z < d/2: 

ay = (z— zi)(z— 3zi+ 2z-y) —{zi — z^'f , (A4) 

for X2 < X < a, z-^ < z < d/2: 

ay — {z — z-Zf' . (A-5) 

For —a < X < a, z > d/2 (above the strip) one has 

Ay{x, z) = Ay{x, d/2) - {Jc/2)F{x, h){z - d/2). (A6) 

The appropriate expressions below the core foUows from 
the symmetry relationship Ay(x^—z) — -~Ay{—x^z). 
Outside the core width one has for x > a and all z in- 
side or close to the strip [for x < —a use Ay{x,z) = 
-Ay{-x,-z)]: 

Ay{x,z) ^ Ay{a,0) + f H^{x' ,Q)dx' -Haxz 

J a 

+ ^g{x,z) (A7) 

where Ay{a,Q) = dH^,-^/2Jc, g = z^/d inside and 
g = |z| — d/4 outside the strip for x < w, g = for 
X > w, and Hz{x,0) from Ref. ^3- Using formulas 
arctanh(l/u) — arctanh(M) + z7r/2 (at \u\ < 1) and 
arctanh(TO) ~ iarctan(u), we may write one single ex- 
pression valid for all — cxd < x < 00 (Re means real part): 



, , Jc { l-a'^/x'^ 

Hzix, 0) = — Re< arctanhW — 

TT V 1 — a 



(A8) 



APPENDIX B: MAGNETIC MOMENT 

From the general definition of the magnetic mo- 
ment of the strip per unit length, (M^jM^) — 
J dx ^ dz{—z,x)Jy{x, z), one obtains the following sat- 
uration values in the two limiting cases: For Ha > — 
(Jc/t^) ln{2ew/d)^ along z, one has j{x, z) = — jc sign(x) 
and = Mf'*' = -jcdw'^ = -Jcw'^. For Ha > Jc/2 
along X, one has j{x,z) = jc sign(z) and = M^^' = 
—jcd'^w/2 — ~Jcdw/2. The reduced magnetic moment 
niz = Mz/Ml'^^ along z depends only on Haz = Ha cos6, 
and is given by Eq. H23I) . The magnetic moment along x, 
Eq. in general depends on both components and has 
to be computed from the current fronts. Explicit formu- 
las for nix = Mx/M^^^ can be obtain for any Ha, but here 
we present them only in the case Ha > Hp = Jc/2 sin 6'. 
Namely, for scenario 1 the formulas of Sec. II yield for 
h > h„ 



pxi 




rrir ~ dx 




Jo 





dx 



1. 
2 

—a 

TT 



1-1 — (h — h^,) shi9 H — arctan 

TT TT 



xVl - 


a2\ 




x^J 


xy/l - 


a2 


Va^- 


X'^ 




— x^ 



COS0 



(Bl) 



with h^, = arcosh(sin0/x)/ cos 9, Eq. (9). For scenario 2 
the formulas of Sec. Ill yield for h> hp 



rux 



r 




/ dx 




'0 





'2 xVT 
— arctan 

TT 



(B2) 



which does not depend on Ha 
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